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^The  i 


i objective  of  this  project  was  to  provide  AP-1 20B  array  processors  with  a program- 
controlled  source  of  15-bit  random-number  words  at  the  rate  of  three  per  microsecond.  A simple 
TTL  circuit  was  implemented  to  do  this.  The  implementation  and  testing  of  a fast  algorithm  to 
generate  uniformly  distributed  random  numbers  is  described.  The  same  shift-register  method  can  be 
expanded  to  include  a fast  RAM  memory  that  will  provide  random  numbers  of  any  given  distribution 
without  slowing  the  controlling  array-processor  program— 
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OBJECTIVE 


Provide  the  AP-1 20B  array  processors  with  a program-controlled  source  of  1 5-bit 
random-number  words  at  the  rate  of  three  per  microsecond. 


RESULTS 

1 . Installation  of  a random-number  generator  on  the  back  panel  of  each  of  five 
AP-1 20B  array  processor  computers  has  been  successfully  completed.  Installation  on  the 
four  remaining  array  processors  will  be  completed  after  acceptance  of  the  hardware. 

2.  A fast  source  of  n-bit  pseudorandom-number  words  can  be  implemented  with  n 
shift  registers  operating  in  parallel,  each  register  supplying  one  bit  to  the  random  number 
word. 

3.  The  system  installed  generates  numbers  with  uniform  distribution.  A memory 
look-up  table  can  be  added  for  generating  numbers  of  any  distribution,  without  loss  of 
speed. 

4.  The  circuit  can  be  implemented  with  standard  TTL  integrated  circuits. 

5.  The  programs  listed  in  this  report,  used  to  debug  the  circuits  (COMPARE2)  and 
to  check  for  digital  noise  (COMPARE  1 and  TIMTEST),  showed  that  the  system  was  suffi- 
ciently noisefree  and  operationally  consistent.  In  one  check,  the  five  generators  were  run 
continuously,  at  high  speed,  for  1 hour.  Each  generated  the  same  number  of  random-number 
words  (about  1010),  and  all  stopped  on  the  same  words. 


RECOMMENDATION 


To  make  the  random-number  generator  more  versatile,  use  the  easily  accessed  word- 
count  register  in  the  array  processor  for  the  following  purposes: 

1 . To  load  the  initial  conditions  into  the  1 5 shift  registers.  This  would  allow  different 
sequences  of  the  given  distribution. 

2.  To  load  a RAM  memory  look-up  table  with  numbers  of  a given  distribution.  This 
would  allow  sequences  of  any  distribution  to  be  generated. 
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INTRODUCTION 


A hardware  implementation  of  the  shift-register  pseudorandom-number  generator 
described  in  NOSC  TN  345*  was  attached  to  the  AP-120B  computer  array  processors  (AP) 
in  the  real-time  torpedo  simulation  system.  Installation  of  a random-number  generator  on 
the  back  panel  of  each  of  five  AP  computers  has  been  successfully  completed.  Installation 
on  the  four  remaining  array  processors  will  be  completed  after  acceptance  of  the  hardware. 
This  document  describes  the  circuit  design  of  the  number  generator,  its  computer  interface, 
the  diagnostic  software  written  to  test  it,  and  the  expansion  of  this  generator  to  a system  that 
generates  numbers  of  an  arbitrary  distribution. 


CIRCUIT  THEORY 

The  circuit  comprises  15  shift  registers  of  16  bits  each.  AH  15  registers  shift  simul- 
taneously (in  parallel)  to  generate  a new  random-number  word.  Bit  0 of  each  register  holds 
the  exclusive  OR  of  bits  7 and  15.  The  15-bit  random-number  word  is  composed  of  bit  1 of 
each  of  the  1 5 shift  registers  (figure  1 ).  A new  number  is  generated  with  only  one  shift; 
whereas  if  but  a single  register  were  used,  about  10  shifts  would  be  required  to  insure  that 
the  numbers  were  not  correlated  with  each  other,  increasing  the  generation  time  by  about  a 
factor  of  1 0. 
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Figure  1.  Random-number  generation  from  15  parallel  shift  registers. 


Each  register  word,  bits  1 thru  1 5,  repeats  only  every  2^  shifts  (a  maximal-length 
sequence).  It  was  observed  that  the  random-number  word  sequence  generated  by  the  1 5 
parallel  registers  is  also  of  maximal  length  when  the  initial  conditions  (noise  words)  given  in 

* NOSC  TN  345,  A Fast  Pseudo-Random  Number  Generator,  by  K Lawrence,  1 February  1978.  NOSC  TNs 
are  informal  documents  intended  chiefly  for  internal  use 
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TN  345  are  used.  These  15  initial  noise  words  were  generated  from  a single  resister  by  load- 
ing bits  1 through  1 5 with  the  sequence  62732  (octal)  and  forming  another  sequence  by 
sampling  bit  1 every  13  iterations.  The  initial  noise  words,  listed  in  figure  2,  are  successive 
1 5-bit  bytes  of  this  sequence.* 
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Figure  2.  Initial  contents  of  the  15  shift  registers. 


SYSTEM  DESIGN 

Each  of  the  nine  AP-120B  array  processors  in  the  simulation  system  will  have  its  own 
noise  generator  attached.  Some  design  criteria  are  as  follows: 

1 . The  AP  request  for  a number  from  the  generator  is  signaled  by  a transition  of  the 
FLAG  3 signal  in  the  AP.  This  flag  is  set  and  cleared  by  the  AP  program  instructions  SFL3 
and  CFL3.  To  change  FLAG  3 and  read  a number  into  the  AP  requires  two  instructions  of 
167  ns  each.  The  generator  therefore  must  be  able  to  produce  numbers  in  about  200  ns  to 
allow  for  transmission  and  AP  memory  setup  time. 

2.  The  FLAG  3 signal  can  be  used  together  with  the  AP  RUN  signal  to  control 
resetting  of  the  number  sequence.  There  are  two  modes  of  reset: 

a.  The  generator  resets  when  the  RUN  signal  goes  low  - ie  when  the  AP 
starts  running  - independently  of  FLAG  3:  RESET  = RUN. 


* From  conversation  wilh  K Lawrence,  NOSC 
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b.  The  generator  resets  when  RUN  goes  low  if  FLAG  3 is  initially  set: 

RESET  = RUN  • FLAG  3.  With  this  mode,  the  number  generator  can  be 
programmed -to  start  on  any  of  the  32767  numbers  in  the  sequence. 

The  desired  reset  mode  is  selected  by  a switch  on  the  AP.  (Switch  2 open  = mode  a:  closed  = 
mode  b.) 

3.  The  port  into  the  AP  is  the  16-bit  HMA2  register.  The  LSB  of  HMA2  is  always 
0,  leaving  the  1 5 MSB  for  the  random-number  bits.  The  noise  generator  is  the  only  source  of 
data  at  this  port.  (AP  device  address  =11.) 

CIRCUIT  DESIGN 

The  block  diagram,  figure  3,  shows  the  interface  with  the  array  processor.  The  shift- 
register  circuit  (figure  4)  uses  the  8-bit  parallel-load  TTL  chip  74198.  Two  chips  are  used  to 
make  bits  1 through  1 5 of  the  register  word.  Register  bit  0 is  the  output  of  the  exclusive  OR 
chip  7486.  There  are  1 5 register  circuits,  each  supplying  one  bit  to  the  1 5-bit  random-number 
word. 


AP-120B  ||  RANDOM- 


NUMBER 


Figure  3.  Random-number  generator  - uniform  distribution  and  fixed  initial  conditions. 


Two  inputs  to  each  of  the  30  shift-register  chips  control  shifting,  (S),  and  loading  or 
resetting,  (L).  Resetting,  by  parallel  loading  of  the  hard-wired  values  at  the  register  inputs, 
occurs  if  the  L signal  is  high  when  S rises. 

The  timing  and  pulse  shaping  circuit  (figure  5)  uses  74121  one-shots  for  timing.  T2 
delays  the  clock  pulse  with  respect  to  the  rise  of  the  load  pulse,  to  satisfy  the  shift-register 
setup  requirement.  The  timing  is  shown  in  figure  6.  The  RUN  signal  is  low  for  the  duration 
of  the  AP  program  and  FLAG  3 is  constant  during  transitions  of  RUN. 


Figure  Shift  register  circuit  (1  of  15). 


CAPACITANCES  ARE  IN  PICOFARADS 

SW2  OPEN:  RESET  - RUN 

SW2  CLOSED:  RESET  = RUN  • FLAG  3 


Figure  5.  Timing  and  pulse  shape  circuit. 


6 


APRUN 

AP  FLAG  3 


1 


■4  b 


:^=x 


M (LOAD)  ' 
(a)  (DELAY)  • 


RESET  = EUR  • j(SW2  OPEN)  + [(SW2  CLOSED)  • FLAG  3$ 

T 


(b)  (SHIFT)  ■ 

© + ® 

(s)  (SHIFT) 

©+® 

© + ® 

@ (SHIFT)  • 


i(LOW) 


(HIGH) 


J 


200 


i nr 


-I  b 
■4  b 


j 150  L 


■4  b 


J 150  1 


-I  t- 


RE5ET  >=  (SW2  CLOSED)  • FLAG  3 


■ 4 b 


H b 


PULSE  LENGTHS  ARE  IN  NANOSECONDS 
Figure  6.  Timing. 


200 


200 


200 


200 


Th 

Tw 


1. 

1-. 


DESCRIPTION  OF  TEST  PROGRAMS 

Three  FORTRAN  programs  were  written  to  test  the  random-number  hardware.  The 
programs  are  entered  into  the  host  computer  (UNIVAC  1 1 10);  these  call  FORTRAN  sub- 
routines which  in  turn  enter  machine  code  programs  into  the  AP.  The  main  programs  are  as 
follows: 

I.  COMPARE2  (compare  to),  which  compares  an  arbitrary  length  sequence  of  AP 
(hardware ^generated  numbers  to  software-generated  numbers  from  the  subroutine  UFORM2. 

II.  TIMTEST  (time-test),  which  starts  the  number  generator  in  a continuous  run  for 
a given  number  of  seconds,  minutes,  or  hours.  Several  APs  can  run  simultaneously  with  this 
program;  at  the  end  of  the  run,  the  APs  can  be  checked  against  each  other  to  see  if  the  final 
random  number  word  in  each  is  the  same. 

III.  COMPARE1,  which  runs  the  generator  over  the  same  sequence  a given  number 
of  times  to  check  for  repeatability. 

To  generate  a continuous  sequence  of  arbitrary  length  in  COMPARE2,  SW2  must  be 
closed  so  that  the  generator  will  not  reset  on  subsequent  calls  to  the  AP  from  the  host  compu- 
ter. For  COMPARE  1 , the  switch  must  be  open  so  that  the  generator  does  reset.  TIMRUN 
will  operate  with  either  switch  setting. 
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COMPARE  1 and  COMPARH2  also  can  list  random  numbers  and  produce  diagnostics 
when  errors  in  the  hardware  numbers  are  detected.  COMPARE2  can  be  used  to  determine 
which  of  the  15  shift  register  circuits  is  faulty,  since  it  outputs  the  correct  word  and  the 
erroneous  hardware  word  in  octal  format.  Written  in  base  2,  the  individual  bits  in  disagree- 
ment represent  faulty  register  circuits.  COMPARE  1 is  conveniently  used  for  listing  since 
SW2  is  normally  ir.  the  open  (RESET  “very  call)  position.  Other  program  options  and  inputs 
can  be  found  under  their  headings,  listed  in  the  appendix.  The  programs  listed  in  this  report 
were  used  to  debug  the  circuits  (COMPARE2)  and  to  check  for  digital  noise  (COMPARE  1 
and  TIMTEST).  In  one  check,  the  five  generators  were  run  continuously,  at  high  speed,  for 
1 hour,  each  generating  the  same  number  of  random-number  words  (about  10 1U);  all  genera- 
tors stopped  on  the  same  word. 

GENERATING  NONUNIFORM  DISTRIBUTIONS 

Random  numbers  of  a given  distribution  other  than  uniform  can  be  obtained  with  a 
table  memory  (look-up  table)  by  using  the  uniformly  distributed  numbers  as  addresses  of  the 
numbers  with  the  given  distribution.  This  is  a direct  implementation  of  the  definition  of  a 
random  variable.  The  equally  likely  outcomes  of  the  experiment  (domain  of  the  random 
variable)  are  the  addresses,  and  the  values  of  the  random  variable  are  the  contents  ot  memory 
The  probability  density  function  of  the  random  variable  - or  its  integral,  the  distribution 
function  — is  determined  by  which  numbers  are  loaded  into  memory. 

The  access  time  of  bipolar  memories  is  about  40  ns  - a time  short  enough  that  a 
memory  look-up  could  be  included  with  the  present  generator  without  slowing  the  random- 
number  generation.  The  maximum  output  rate  will  still  be  one  random-number  word  every 
two  instruction  cycles  (334  ns).  The  memory  can  be  PROM  or  RAM.  If  RAM  is  used,  differ- 
ent distributions  are  programatically  available.  A convenient  register  is  included  on  the  AP 
237  board  through  which  the  RAM  can  be  loaded:  the  word  count  (WC)  register,  adjacent 
to  the  HMA2  register.  The  shift-register  initial  words  can  also  be  loaded  via  the  WC  register, 
to  make  available  different  sequences  of  uniformly  distributed  numbers. 

The  numbers  to  be  loaded  originate  in  the  1110  host  computer.  They  are  read  into 
the  AP  memory  (MD),  then  passed  one  at  a time  through  the  WC  register  to  the  RAM  or  to 
the  shift-register  inputs.  Figure  7 is  a block  diagram  of  a design  that  allows  the  loading  of 
both  a RAM  and  the  initial  register  words  and  provides  two  programmably  selectable  outputs 
- either  the  uniformly  distributed  words  from  the  shift  registers  (the  RAM  addresses)  or  the 
nonuniform  words  from  RAM.  Although  uniform  distributions  are  available  from  RAM. 
sequences  from  RAM  are  limited  by  its  size  (about  4k  words). 

The  WC  register  presents  a control  word  before  each  data  word.  Control  words  are 
distinguished  by  those  which  have  1 as  the  most  significant  bit.  The  other  bits  in  a control 
word  are  used  to  indicate  where  the  following  word  of  data  should  be  put,  to  enable  busses, 
and  to  set  the  output  select  latch.  A data  word  is  loaded  into  a shift  register  or  RAM  after  a 
delay  initiated  by  the  falling  edge  of  WC  bit  0.  Figure  8 shows  the  bit  assignments  of  the  WC 
register. 

The  full  circuit  shown  in  figure  7 requires  three  circuit  boards,  each  about  the  same 
size  as  the  present  one  (6  by  10  inches).  They  are  a shift  register  board,  a RAM  board,  and  a 
control  board. 
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Figure  7.  Random-number  generator  - arbitrary  distribution  and  initial  conditions. 
bit  o 1 


t5 


Figure  8.  Bit  assignment  in  AP  WC  register. 

CONCLUSIONS 

1.  A fast  source  of  n-bit  pseudorandom-number  words  can  be  implemented  with  n 
shift  registers  operating  in  parallel,  each  register  supplying  one  bit  to  the  random  number 
word. 
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2.  The  system  installed  generates  numbers  with  uniform  distribution.  A memory 
look-up  table  can  be  added  for  generating  numbers  of  any  distribution,  without  loss  of  speed. 

3.  The  circuit  can  be  implemented  with  standard  TTL  integrated  circuits. 

4.  The  programs  listed  in  this  report,  used  to  debug  the  circuits  (COMPARE2)  and 

to  check  for  digital  noise  (COMPARE  1 and  TIMTEST),  showed  that  the  system  was  suffi- 
ciently noisefree  and  operationally  consistent.  In  one  check,  the  five  generators  were  run 
continuously,  at  high  speed,  for  1 hour.  Each  generated  the  same  number  of  random-number 
words  (about  1 0 ),  and  all  stopped  on  the  same  words. 


RECOMMENDATION 

To  make  the  random-number  generator  more  versatile,  use  the  easily  accessed  word- 
count  register  in  the  array  processor  for  the  following  purposes: 

1.  To  load  the  initial  conditions  into  the  1 5 shift  registers.  This  would  allow  different 
sequences  of  the  given  distribution. 

2.  To  load  a RAM  memory  look-up  table  with  numbers  of  a given  distribution.  This 
would  allow  sequences  of  any  distribution  to  be  generated. 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

tl 

12 

13 

14 

15 
14 

17 

18 
19 
<0 
21 
22 

23 

24 

25 
24 

27 

28 

29 

30 

31 

32 

33 

34 

35 
34 

37 

38 

39 

40 

41  2 

42 

43 

44  ( 

45 
44 

47  C 

l 

49 

50 

51 

52 

53 
5* 


PROGRAM  T 1 7 IE  ST 

TO  TEST  THE  RAND  OH  NUMBER  6ENERAT0R  ATTACHED  TO  THE  AP120B 
ARRAY  PROCESSOR.  RUNS  THE  NUMBER  GENERATOR  IN  NESTED  LOOPS, 

EACH  LOOP  READING  2 BOROS  INTO  THE  AP  MEMORY  (MD  Q)  . TO  REAR 
ONE  BORD  REQUIRES  2 I NSTRUC TI ONS *3 34  NS.  THE  TOTAL  RUN  TIME  IS 
CLOSE  TO  THE  PRODUCT  OP  THE  THREE  LOOP  COUNTERS  (N215,  N230,  N2451 
AND  334  NS. 

FOR  SHORT  RLNS  (LESS  THAN  1 SEC.)  THE  AP  BILL  HRITE  THE  FINAL  UORi  TO 
THE  1110  HOST  COMPUTER.  FOR  LONGER  RUNS  THIS  PROGRAM  MUST  8E  TERMINATED 
AND  THE  FINAL  WORD  READ  FROM  THE  AF  (FROM  HD  0) 

WITH  THE  SYSTEM  PROGRAM  BUG.1  OR  BIG. 2.  THE  TERMINAL  INPUTS  ARE: 

9X9 T BUG.1 

(ARRAY  PROCESSOR  NO.,  0-9,  11) 

E (EXAMINE) 

MO  (MEMORY) 

0 (ADDRESS  ZERO) 

X (TO  TERMINATE  P R 06 RAM ) 


INPUTS 

N215-  PRIMARY  LOOP  COUNT  (327(6  MAX) 
N230«  2-ARY  " (32767  MAX) 

N2  *5=  3-ARY  " 

TO  RUN  1 SEC  SET  N21*=32766 
N23C=92 
N24  5 =1 


TO  RUN  1 PIN  SET  N215*32766 
N230 - 92 

N24  5*60 


TO  RUN  1 HOUR  SET  N215*32766 
N230*92 
N24  5 =3600 

ONE  HOUR  REPRESENTS  32766*92*3600*  10**10  SHIFTS  OF  THE  NOISE  REGISTERS. 
BILL  LAFOND  APR  78 

FORMAT ( 15 ) 

FORHATOS  ,06) 

CALL  APCLR 

REAMS, 1)  N21 5 ,N230,N2 45 

SET  FLA6  3 

CALL  TIMRUN  (2  ,0,0) 

CALL  APHR 

START  TIMED  RUN 

CALL  TIHRUN(N215,N230,N24S) 

CALL  APHR  

CALL  APGET (NFINAL, 0,1,0) 

CALL  APBD 

HR 1TE  (6 ,2)  NFINAL  

STOP 

END 


1.  ST1TLE  HARDWARE  RANDOM  NU*HtR  TIME  TEST  PRC6R AM  1 1 MR UN 

2.  SINT  NT  TIMRUN, 3 

3.  ■ runs  RANDOM  NO.  GEN  IN  NESTED  LOOPS  SO  THAT  IT  HILL  RUN  A PRFDEIRMINED 

N . " LENGTH  OF  TIME,  E .».  , HOURS,  TO  TEST  IF  ALL  OF  THE  ARRAY  PROCESSORS 

5.  " ARE  RUNNING  CONSISTANTLV  AND  THE  SAME. 

6.  - SFADO*  COUNt  UF  Tc  <**15  - 2 ITERATIONS  * 11  rS.  PAX. 

7.  • SPAM*  ■ 4.»»30  - 1 " - < MINUTES  I»  SPU  AND  SP1*HAX 


12 
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F.  " SFAC2*  " <**45  - 1 " * 1 NCOS  IF  *10  ANO  OT  HERS  Ml. 

V.  - TOTAL  TIME*  TNF  PRODUCT  J 34  NS  . * ( SP AO  C ) * ($ p A D1  ) • CSP A D?  ) 


lb. 

m 

BILL  LAFOND  APR  It 

11. 

m 

12. 

1IMRUN:  NCP 

13. 

DB* 

l ERO;  LDAPS 

14. 

D6*0;  LOSPI  15 

15. 

MOV 

15,15;  SETMA;  C6*ItR0;  M 1 <D B 

1e. 

DEC 

0;  INCMA 

17. 

MOV 

0,3 

lb# 

MOV 

1,4 

IS. 

LDDA;  DB*  1 1 

20. 

LL: 

IN;  D6*INBS;  MKOB;  aRTLPN;  DECMa 

21. 

DEC 

C;  CF  L3;  INCMA 

22. 

IN; 

DB*InBS;  MKDB;  wRTLHN;  DECMA 

23. 

DEC 

0;  B6T  LL;  SFL3;  INCMA 

24. 

DEC 

1 

25. 

BG1 

LL;  NOV  3,0 

26. 

DEC 

2 

27. 

BCT 

LL;  MOV  4,1 

RETURN 

29. 

SEND 

1 

2 

3 

A 

5 

6 
7 
• 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 
28 

27 

28 
29 
3C 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


PROGRAM  C OP  FA  RE  1 

TO  TEST  THE  A A NO OH  NUHBE R GENERATOR  ATTACHE*  TO  THE  AP  1206 
ARRAT  PROCESSOR.  REAPS  NS  SEQUENCES  Of  RANDOM  NO.'S.  EACH 
SEQUENCE  NR  HOADS  10N4.  SWITCH  2 ON  THE  RANDOM  NO.  6ENERAT0R  SHOULD 
BE  OPEN  SO  THAT  THE  NUMBER  Oi HI  BATCH  RESETS  ON  EVERT  CALL  TO  THE  AM. 
THE  PROGRAM 

1.  COMPARES  EACH  SET  TO  THE  PROCEE  DIN6  BV  C0HPAR1N6  TEST  SUNS 
COMPUTED  IN  THE  AP.  IF  THE  SUMS  ARE  DIFFERENT  TNI  HORDS  AH  READ 
FROM  THE  AP. 

2.  CHECKS  THAT  SUCCESSIVE  WORDS  ARE  DIFFERENT. 


INPUTS 

NS *N0 . SETS  TO  BE  COMPARED  (L.E . 999°9> 

N8*N0.  RANDOM  NO.'S  IN  EACH  SET  <L.t.  20001 

ND I A6 '1  TC  LIST  RANDOM  NO.'S  AND  DIAGNOSTICS.  TO  LIST  RANDOM  NO.'S 
SET  N S* 1 , NDIA6*1,  NR*LEN6TH  OF  SEQUENCE  DESIRED 
NP A SS * 1 TO  MARK  END  OF  EACH  COMPARISON 


BILL  LAFONO  APR  76 
.INTEGER  R (2  COli  > . w(20CC) 

FORMAT  CIS) 

FORMATC10C5l.06>> 

F0RMAT(/T5,'SET  NO.  1') 

FORMATES, 'DIFFERENCE  FOUND  IN  RANDOM  NO.  '.IS, 

C ' OF  SET  NO.  ',15 .) 

FORMATCTS.'PASS  ',15,'  FINISHED'/! 

F0RMATCT8, 'PRESENT  WORD*  ',06,'  LAST  WORD*  ',06,1 
F0RMATCT5, '«•••**  aORO  ',16,'  IS  SAME  AS  LAST  ONE  - ',06, 
C 'IN  PASS  ',It,I 

FORMA ? C T5 , * 16  LSB  OF  SUM  * ',06) 

0 F0RMATC/T5 ,'END  TEST') 

CALL  APCLR 


I 

1 


1 
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*2 

RE  AO  ( 5 « 1 ) NS,NR,NDIA6,NPASS 

43 

C 

44 

(ALL  HANOI (C,1,NR,4094> 

45 

CALL  APwR 

44 

CALL  AP6E T (ISOM) ,4094 .1.0) 

47 

CALL  APWO 

48 

If  (NO  1A6.M  .0)  w R I TE  ( 4 ,8  ) LSUM1 

49 

c 

SO 

c 

MITE  fIRST  SET  Of  RANDOM  NO.'S  IN  OCTAL 

51 

t 

52 

CALL  AP6iT(fi,0,NR,C) 

53 

CALL  APWO 

54 

J1=1 

55 

DO  15  1=2, NR 

54 

If  (R(I).EG.R(l-D)  •R1TE(4,7)  I,R(1),J1 

57 

15 

CONTINUE 

58 

If  END  IA6.NE  .0)  bR1TE(4,3) 

59 

If (NO 1AG.NE  .0)  WRI TE  <4  ,2  > ( R ( 1 ) , I = 1 , NR ) 

40 

If  (NPASS.NE.O)  W R1 TE (4  ,5 ) J1 

41 

c 

42 

If  (NS  .E9.1 ) 60  TO  999 

43 

DO  100  Kl=2,NS 

44 

CALL  HANOI (C,1, NR, 4094) 

45 

CALL  APWR 

44 

CALL  APGE T (LS0H2.4C94 ,1,0) 

47 

CALL  APWO 

48 

If  (NO 1A6.NE  .0)  wNl TE (4  ,8 ) L S UH 2 

49 

LD I f = LSUM 1 - LSUM2 

70 

If(LDIf.EQ.C)  60  TO  99 

71 

c 

RETRIEVE  PRESENT  S')  THAT  DOESN'T  AfaREE  WITH  LAST  ONE 

72 

c 

73 

CALL  AP6ET  (W.CtNR.O) 

74 

CALL  APWO 

75 

If  (W(I).NE  .R(1))  WRITE (4,4)  J 1 ,K 1 

74 

If  (W(I).NE.R(D)  WRITE  (4, 4)  W(1),R(1) 

77 

DO  20  1=2, NR 

78 

If  (W(l).ER  .W(i-1 ))  WRITE  (4,7)  I,W(1),K1 

79 

If  (W(l).EQ.W(l-D)  60  TO  21 

40 

If  (M(I).E«>R<I)>  60  To  20 

81 

WR  IT  E (6  ,4  ) 1,K1 

82 

WR ITE (4  ,4  ) W(1),R(1) 

83 

60  TO  21 

84 

20 

CONTINUE 

85 

21 

If (NDIA6.NE  .0)  WRITE  (4,2)  (W ( I ) , I = 1 ,NR ) 

44 

DO  22  1=1 , NR 

87 

22 

R(  I ) = W ( I ) 

88 

99 

If  (NPaSS.NE.O)  WR | T| (4  *5 ) K 1 

89 

LS  UM1 =L  SUM  2 

50 

100 

CONTINUE 

VI 

999 

CONTINUE 

92 

WRITE  (4,10) 

93 

STOP 

94 

END 

1. 

STITLE  HARDWARE  RANDOM  NUMBER  TEST  PROGRAM  RAN 01 

2. 

SENTRY  RAND  1 »4 

3. 

S EXT  VCIR 

4. 

m 

BILL  LAfONO  APRIL  7e 

5. 

m 

fOR  NO  SC  SPECIAL  MOD IFICATION  TO  ATTACH  HARDWARE  RAnOO*  NO* 

4. 

m 

GENERATOR  TO  THl  AP  1200. 

7. 

m 

LOADS  THE  FIRST  (SPAD2 ) LOCATIONS  IN  HD  MEMOR T WITH  RANDOM  NO.'S 

8. 

m 

AND  CALLS  SUBROUTINE  SVE1  TO  SUM  THE  NC.'S  ANO  FLACES  THE  SUM 

9. 

m 

IN  LOCATION  (SPADj)  IN  MD . 

14 
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<c. 

a RANH  IS  CALLED  BY  PROGRAM 

COMPAREl. 

11. 

■ switch  z on  the  random  no. 

gen 

. SHOULD 

12. 

OR 

13. 

" SPADO:  fa A SE  ADDRESS 

14. 

" SFADl*  INCREMENT 

15. 

" SPA02*  no.  of  random  no.'s 

It. 

■ SPAD3-  MEMORY  LOCATION  OF 

TEST 

SUM 

17. 

R AND  1 : NOP 

IE. 

MOV  0,7 

19. 

MOV  1,10 

2C. 

MOV  2,11 

21. 

MOV  3,12 

22. 

OB=ZERU;  LDAPS 

23. 

mov  0,e 

24. 

J SR  VC  LR 

2S. 

MOV  6,0 

2D. 

CFL3  "INITIALIZE  FLAG  3 

27. 

NOP 

28. 

LDDA;  DB= 11 

29. 

L DSP  1 6;  0B=177777 

30. 

MOV  D,t;  SETHA 

31. 

MOV  2,0  "SPt  IS  COUNTER  REG 

32. 

DEC  6 

33. 

LL:  IN,  DB  = INBS;  PKOB;  1NCMA 

; WRTLMN 

34. 

S Fl3;  DEC  D 

35. 

IN;  DB-1NBS;  MKDB;  INCMA;  WRTLMN 

3D. 

CFl3;  DEC  D;  86T  LL 

37. 

NOP 

38. 

NOP 

39. 

JSR  SVtl 

40. 

MOV  7,0 

41. 

MOV  10,1 

42. 

MOV  11,2 

43. 

MOV  12,3 

44. 

RETURN 

45. 

SVEI:  NOP  "ADOS  CONTENTS  OF 

HD  - 

SPO=  BA! 

44. 

NOV  0,0;  SETHA  " .SP2=N0. 

OF  WORDS,  SP 

47. 

LDSP1  D;  DB=0 

48. 

MOV  2,2;  Db'ZERO;  DPYCOXDB 

"FLOAT.  0 TO 

49. 

LOOP:  LDSPI  5;  QB  *HD ; UEQ  LOOPEND 

SO. 

ADD  1,0;  SETHA 

51. 

ADD  5,0 

52. 

DEC  2;  BR  LOOP 

S3. 

LOOPEND:  NOP 

54. 

mov  0,0;  db  = spfn;  opaCGXDB 

SS. 

DB*DPIC0);  DPYCOXDB;  WRTLMN 

SC. 

MOV  3,3;  SETHA;  MXDPy<0) 

SI. 

RETURN 

S4. 

SEND 

INCREMENT 


1 

2 

3 

4 

5 

6 
7 
6 
9 

10 

11 

1? 

13 

14 


PROGRAM  CL'MNti 

TO  T I ST  THE  RANDOM  NUMBER  GENERATOR  ATTACHED  TO  THE  AP  1208 
ARRAY  PROCESSOR.  THE  PROGRAM  READS  A CONTINUOUS  SEQUENCE  Of 
RANDOM  NUMBERS,  ESSENTIALLY  UNLIMITED  IN  LENGTH,  BY  READINC  NS 
SETS  OE  NR  VORDJ  CRANDOM  NO.'S)  E*CH. 

41 TH  THIS  SEQUENCE  , THE  PROGRAM 

1.  SEARCHES  EOM  THE  REPETITION  Cf  A WORD  Of  6IVEN  SEQUENCE  NO.  IN  SET 
1,  AND  UUTPUTS  THE  WORD  AND  THE  SEQUENCE  NO.  OF  THE  REPEATED  WORD. 
ANY  WORD  REPEATS  EVERY  2«*15  ITERATIONS. 

2.  SEARCHES  FOR  A WORD  AND  OUTPUTS  ITS  SEQUENCE  NO. 


15 


VIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
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rW 


|5  C 3.  COMPUTES  THt  RANDOM  NO.'S  WITH  THE  SUBROUTINE  UFORH2  AND 

U C COMPARES  THE  CONFUTED  SEQUENCE  TO  THE  HARDWARE  (AP)-OENERATED 

17  C WORDS. 

16  C 

19  C TO  COMPARE  THE  SAME  SEOUENCE  .1IH  ITSEIE  N TIMES  SEE  THE  P R CORAH 

20  C COM  PA  RE  1 • BILL  LA  f ON  D APR  7 2 

II  C 

22  C INPUTS 

23  C 

24  C NS=NO.  SETS  TU  d E bENERATEB  (L.E.  VV999) 

25  C NR  = NO • RAkDOM  NO.'S  IN  EACH  JET  (L.E.  2000) 

26  C NO  I At  = 1 TO  LIST  WORDS.  Ci  TO  NOT  LIST 

27  C NUDlAb  =1  TO  PRINT  bCROS  FROM  SUL.  UF0RA2 

28  C NPASS  * 1 TU  MARK  ENOS  OF  LOOPS 

29  C NCHk  = SER.  NO.  OF  kORD  IN  FIRST  SET  TO  BE  COMPARED  TO  ALL  OTHERS 

30  C =G  TC  SNIP 

31  C NCMPR  =1  TC  COMPARE  UFO RH 2 (SOFT»»nE>  WORDS  10  HARDNARE  NORDS 

32  C =0  TO  SMP 

33  C NF I NO  = (OCTAL)  «OND  TO  BE  FOUND  - OUTPUTS  SEO.  ANO  PASS  NO. 

34  C = 0 TO  SKIP 

35  C 

36  INTEGER  R (2C00l.Uk  (2000 

37  C 

38  1 FCRMATU5) 

39  t FORMAT (10(5*,  06 >) 

40  3 F0RMAT(/T5,'SET  NO.  1') 

41  4 F0RNAT(/T5  , 'RANDOM  NUMBERS  - SET  ',15.) 

42  5 F0RMAT(T5 , 'PASS  ', 1 5 , ' FINISHED'/) 

43  7 FORMAl(T5,'bORD  ',It,'  IS  SAME  AS  WORD  '.16,'  * '06) 

44  8 FORMAICTS,'??????  SOFTbARE  WORD*  ',06,'  HARDWARE  WORD*  '. 

45  C 06,'  bORC  NO.  ',16) 

46  1C  FORHAMTS.'ENW  SOFT  - HARD  COMPARE') 

47  11  FORMATC06) 

48  12  F0RMAT(T5,'SERUENCE  NO.  ',16.'  IS  RANOOM  NO.  '.06) 

49  C 

50  CALL  APCLR 

51  RLADC3.1)  NS.NH.MOIAG.NUOIAG.MPASS  ,NCHK .NCMPR 

52  RE  AD  (5, 11  > NF1N0 

53  C 

54  CALL  RAND2(C, 1,2.1) 

55  CALL  AP.R 

56  CALL  RAND2  (C.I.NR.O) 

57  CALL  APWR 

58  NR  1 = N RM 1 

59  CALL  AP6E  T (R.G.NR1.0) 

60  CALL  APwD 

61  C 

62  C WRITE  FIRST  SET  OF  RANDOM  NO.'S  IK  OCTAL 

63  t 

64  J1=1 

c5  IF  (NCHK.EO.C)  GO  TC  19 

66  CkbOrR(NCHK) 

67  DO  15  1-1  , RR 1 

68  IF  ((CKWD.h£.R(l)>.OR.(NCHK.E«.l)>  GO  TO  15 

69  URI1E(et7)  l.hCHK.R(l) 

7G  15  CONTINUE 

71  19  CONTINUE 

72  IF  (NO  IA6.NE  .0)  aRlTE(6,4)  J 1 

73  IF (NO IA6.NE  .0)  »R1 TE (6  ,2 ) (fl  (1 ) .1 *1 .NRl ) 

74  |F (hCHPR.EO.C)  bO  TO  3C 

75  CALL  UF0RH2 (UR. NRl ,1 .NUDIA6) 

76  DC  32  I *1  . ARl 

77  IF  (UR(I).NE.R(D)  bkITE(6.8)  UR(I).R(D.I 

78  32  CONTINUE 

79  IF (NPASS. NE  .0)  «R1TE(C.1G> 

80  30  CONTINUE 

el  IF  CNF 1ND.ER.0)  bO  TO  40 


■ 


IBIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
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b2 

00  b 2 1=1,  NR1 

13 

If  (S(  D.Eo.NFlND)  bRI  TE<4,  12)  I.NP1ND 

eb 

* 2 

CONTINUE 

05 

bC 

CONTINUE 

Ob 

IMNPASS.NE.O)  • R 1 TE  ( t ,5  > J1 

a7 

C 

68  ' 

If (NS.E0.1)  GO  TO  559 

89 

00  100  Kl =2  ,NS 

>0 

Kb  =MPl«  (Kl-2) «NR 

51 

CALL  SAND2 (C,1,NR,0) 

92 

CALL  APbR 

53 

C 

5b 

CALL  APGET  (fc,U,NR,C) 

55 

CALL  APUD 

98 

IP (NCHK.EB.C)  GO  TO  21 

57 

DO  17  1=1, NS 

58 

NS  EB*NA ♦! 

59 

IF  (CK.O.EB.RU))  UR  I TE  (6,7)  N S EU ,N CHK ,R  ( I ) 

ICO 

17 

CONTINUE 

1C1 

21 

CONTINUE 

1C2 

If  (N01AG.NE  .0)  wRITECC.b)  * 1 

1C3 

IF  (NDIA6.NE.0)  bRlTt(t,2)  (STD, 1 = 1, NR) 

1Gb 

If (ncmpr.eb  .0)  GO  lo  2C 

1C5 

CALL  UFORN2«UR,NR,C,NOOI AG) 

108 

DO  22  1=1, NR 

1c7 

NS  E«  = Kb«I 

ICf 

IF  CUR  (I).NE.RCI))  SR1TE(6,2)  uR<  I ) ,R  (I)  ,NSEO 

1C9 

22 

CONTINUE 

110 

IF (NPASS.NE .0)  ■RITL(6,10) 

111 

20 

CONTINUE 

112 

IF  (NF  1N0.ER  .0)  GO  TO  5 0 

113 

00  52  1*1  , N R 

11b 

IF<*<]),EO.NFJND)  GO  TO  51 

115 

52 

CONTINUE 

118 

GO  TO  SO 

117 

51 

NS£a*Kb«i 

118 

UN ITE (8 ,1 2 ) NSE8 ,NF I NO 

119 

5C 

CONTINUE 

120 

99 

IF (NPASS.NE  .0)  NRITE(«,S>  *1 

121 

100 

CONTINUE 

122 

99V 

CONTINUE 

123 

STOP 

12b 

ENO 

, 


1 

2 

3 

Sot-ttOlTlNE  IFCHF.2  <F  ,NS,NHES,NitIA6> 

i 

( 

GENERATES  LSIFORNLT  0 I S T N1  B UTE  S ft  A NOON  NUMBERS  f RON  A 

b 

C 

15  bIT  SHIM  REGISTER  (SEE  NOSC  IN  3b5  bT  K.  LAURENCE) 

5 

C 

AND  APPENDS  A 0 TO  THE  LSB  OF  THE  NUMBER  UORO  TO  CAUSE 

< 

c 

AGREEMENT  WITH  THE  HARDWARE  GENERATED  UCROS  IN  THE  AP120  fc 

7 

c 

ASR  AT  PROCESSOR  (NMA  BIT  IE  IS  A l » ATS  0). 

8 

9 

c 

c 

OUTPUT 

1C 

c 

KitNMT  OF  RANOOM  NUMffRS  (INTEGER  OCTAL)  DIMENSIONED  IN  CALLING 

11 

c 

12 

c 

PROGRAM. 

13 

( 

1b 

c 

INFUT 

15 

c 

16 

c 

SM  NO  OF  NLMOEhS  TO  BE  GENERATED 

17 

c 

NS  E S = 1 TO  PISET  TP.I  St  BONCE,  0 TO  NOT 

17 

c 

ND1AG=1  TO  PRIM  DIAGNOSTICS 

19 

c 

17 


WIS  MGS  15  BXST  QUALITY  FRACTICABLf 
FROM  COPY  FURNISHED  TO  DDC  


20 

C 

ADAPTED  FRCP  THE  ORIGINAL  BT  A ARE N LA. PENCE  - BILL  LAFOND 

21 

C 

2? 

C 

23 

DIMENSION  INUMCI5)  ,NUM15> 

2A 

INTEGER  >(1) 

25 

C 

SHIFT  SAMPLINb  bY  13  FOR  SEED  NOS. 

26 

DATA  INUM/C21 77w ,0137 i 0,01 072  7,070626  ,054  726, 

27 

c 

041 76 2,04 61.N, 00501 5 ,0207 00, 0141 fO, 

26 

c 

01 3074, oC 52 *2,04 2u* 6 .055276 ,054435/ 

29 

C 

iO 

1 

FURMAT(/T5,  'UNSHIFTED  15  BIT  INITIAL  CON D IT  I ON  W ORD S FROM 

31 

c 

,'UF0RM2'/> 

32 

2 

FURKAK/T5  ,'SGF1VARL  GENERATED  RANOOP  NUMBERS'/) 

33 

3 

FoRMAT(10(5x,C6)> 

34 

C 

35 

IF  CNRES.EW.O  GO  10  IS 

36 

DO  20  1*1,15 

37 

2 C 

NuMII >*1NUP(I) 

38 

15 

continue 

39 

IF ((NDIAG.EG.uI.OR.CNRES.EO.d)  go  10  18 

40 

•R  ITE  (6  ,1  ) 

41 

WRITE  <6,  3)  (NLM  ( I ) , 1 * 1 , 1 5 ) 

42 

L 

43 

C 

44 

18 

DO  ICO  Kl*1,NR 

45 

C 

MULTIPLY  bY  T.O  TO  SHIFT  ONE  blT  CPMA15*0) 

46 

R(tCl)=NUMU>*t 

47 

1 7 *NUP ( 7 ) 

48 

I15=NUMCT5) 

49 

NE  w*X  OR (17,115) 

50 

DO  1C  1=15, 2,-1 

51 

NUM ( I )~NUM ( 1— 1 ) 

52 

10 

CONTINUE 

53 

NUM Cl )*NE  W 

54 

100 

CONTINUE 

55 

IF  (ND1A6.NE.0>  «RiTE(t,2) 

5t 

IF  (ND 1A6.NE .0)  ■R1TEC6.3)  ( A ( I ) , I * 1 , NR  ) 

57 

RE  TURN 

58 

END 

1.  (TITLE  HARDWARE  RANDOM  NUMBER  TEST  PROOR  AM  RA  NO! 

..  (ENTRY  RAND2.4 

3.  (EAT  VtU 

A.  “ BILL  LAFGND  MARCH  7o 

5.  " TOR  NO  SC  SPECIAL  MODIFICATION  TO  ATTACH  HARDWARE  RANDOM  NO. 

6.  " 6CNLRAT0R  TO  THE  AP  WOO. 

7.  “ LOADS  THE  FIRST  (SPAD2 ) LOCATIONS  IN  Mt  MEMORY  WITH  RANDOM  NO.'S 

0.  " IF  FLAG  3 IS  SET  kHE  N THIS  SUBROUTINE  IS  CALLED  THE  SEQUENCE 

V.  * 1$  INITIALIZED.  OTHERW12E  THE  SEQUENCE  CONTINUES  FROM  THE 

1C.  " LAST  mORO  GENERATED.  SWITCH  2 ON  THE  RANDOM  NO.  6 EN.  SHOULD  BE  CLOSED 

11.  • THE  FORTRAN  SUBROUTINE  RAND?  IS  CALLED  BY  PR06RAM  COMPARE?. 

It.  ” SPADO*  t A SE  ADDRESS  IN  MD  MEMORY 

13.  " SHARI*  MEMORY  INCREMENT 

1A.  * SPA12*  NO.  OF  RANDOM  NO.'S  TO  BE  GENERATED 

IS.  “ SPAD3*  RESET  F L *6  CSETS  F L 3)  > 1*RESET  NEAT  CALL,  0*  NOT 

It. 

17.  N AND2 : NOT 

la.  MOV  0,7 

19.  MOV  1,10 

2 C . MOV  2,11 

21.  CB*Z  LRU;  LDAPS 

?2  . M OF  ti|t 


18 


2i. 

INC  2 

24. 

J SR  VCIR 

25. 

DEC  2 

26. 

MOV  t,J 

27. 

HOP 

21. 

LDDA;  68=11 

25. 

l.  DSP  1 t;  Du=1?77 

3C. 

KOV  t,t;  SETNA 

31. 

KOV  2,6  “SPE  IS 

32. 

DEC  t 

33. 

b Ft3  LLl 

34. 

LL:  IN;  6b=lNBS; 

35. 

SFL3;  DEC  E 

3e. 

Lll : IN;  fb  = 1NBS , 

37. 

IFL3;  DEC  6;  B6 

3t  • 

NOP 

3 V. 

NOP 

4C. 

KOV  7,0 

41. 

KOV  10,1 

42. 

NOV  11,2 

43. 

NOV  3,3  "CHECK 

44. 

E»EU  NORES 

45. 

SM.3 

4t. 

NORES:  NO F 

47. 

RETURN 

48. 

SEND 

0*18  PA©  K BEST  QUALITY  PRACTICABLE 
COPY  FURNISHED  XOLDC  


CCUNTEK  RE6 


P KliB;  1NCNA;  HRTLNN 


PI <0b ; 
IL 


INCPA;  URTLHN 


if  n E$t t n t k t Tike 


